
for (yearx in c(1950:1969,2009:2011)) 

	{
	
	t.maxfiles = maxfiles[which(substr(maxfiles,nchar(maxfiles)-14,nchar(maxfiles)-11)==yearx)]
	
	for (tfile in t.maxfiles)
	
		{
	
		tname = basename(tfile) #Get the just the name of the ASCII file from the file list
	
		tname = gsub('.asc.gz','.sh', tname) #Use tname to create the name for the .sh file
	
		arg.tfile = paste('tfile="',tfile,'" ',sep='') # Define the argument for tfile 
	
		zz = file(tname,'w') # This command creates a file with the name 'tname' The 'w' command opens the file and writes the following lines of text into it
		
			cat('##################################\n',file=zz) # First blank line of .sh file ### Note the use of '\n', any text following this will be on the next line
			cat('#!/bin/sh\n',file=zz) # First command line for HPC
			cat('cd $PBS_O_WORKDIR\n',file=zz) # Sets the directory to write out summary files at the above directory 'tmp.pbs'
			cat("R CMD BATCH --no-save --no-restore '--args ",arg.tfile,"' ",script.file,' ',gsub('.sh','',tname),'.Rout \n',sep='',file=zz) # This line tells R not to save to the .Rdata file and also not to restore the environment from the .Rdata file, very important when batching or you will get multiple fuckups, also defines the arguments, says where script to run is located, and defines the name for the R.out file
			cat('##################################\n',file=zz) # Last blank line of .sh file
	
		close(zz) # Closes the .sh file that you've just written
	
		# The follow line submits the job to the HPC, leave it commented out until you're ready to RUN all the jobs\
		# Always run the script first without removing the comment, this will still generate the .sh files for you to check before running
		# Once you're happy with that then submit a single .sh file and make sure it runs through to completion before batching out the rest with the following line
		
		system(paste('qsub -m n -l pmem=120 ',tname,sep=''))
		system(paste('sleep 1'))

		}
		
	}

# Done
